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Abstract 

We consider the cosmological system with two in¬ 
teracting fluids: dark energy and dark matter, in a 
homogeneous and isotropic universe with dissipation. 
The modified gravitational equation for dark matter is 
solved. The analytic representations for the Little Rip, 
the Pseudo Rip, and the bounce cosmology models with 
dissipation are obtained in terms of the thermodynamic 
parameters in the equation of state. We analyze the cor¬ 
rections in the energy density for dark matter, in view 
of the dissipative processes and the coupling with dark 
energy. 

Keywords dark energy, dark matter, viscous cosmol¬ 
ogy 


1 Introduction 


According to the modern cosmological picture the uni¬ 
verse is filled with a negative pressure fluid (dark 
energy) which accounts for about 73% of the total 
mass/energy of the universe and only 27% o f a combina¬ 


tion of dark matter and baryonic matter; cf. iRiess et al 
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^ Afhlation 


( 1998ll : Perlmutter et al. ( 1999ll . One phenomeno¬ 
logical approach to describe the dark energy is to 
assume that this fluid obeys an equation of state 
with constant thermodynamic parameter w = pjp, 
where p is th e dark energy de n sity and p is t h e dar k 
pressure; cf. Li et al. ( 2^011 1: Bamba et al. ( 2012ll : 
Noiiri and Odintsov ( 20071) . It is known from observa¬ 


tional data that in the current universe the value of w 
lies near, o r pro bably below, —1, w = —l.Odlg;)*®; cf. 


Kowalski (120081) . The region where w < — 1 is called 


the phantom region. 

Various cosmological models for the evolution of the 
universe have been presented. They aim at explaining 
the acceleration era of the universe, and they include 


Big R 

ip Caldwell et al. ( 

0 

0 

Noiiri and Odintsov 

(2004, 

2003)1. the softer Little Rip Framoton et al. 

(2011) 

: Brevik et al. (2011): 

Framoton et al. (2012)]. 


the P seudo Rip iFrampton et al. (l2012ah ]. the Quasi 


Rip IWei et al. (12012111 , a n d the bounce cosmol- 


BrandenbergeiL (I2OII ); NoyellcLaifo BerghaMa. 


ogy _ n—i I_n—rn_ 

(2008): Bamba ( 2013 ): ICai et al. ( 2011 ): Brevik and Timoshkin ^ 

( 20151) ]. The Little Rip and the Pseudo Rip phenom¬ 
ena, and the bounce cosmology pictured as an inho¬ 
mogeneous nonviscous dar k fluid coupled with dar k 
matter, were considered bv iBrevik et al. (120131 . 120141) . 
Cosmological models in which the modification of grav- 


Mvrzakul et al. (2013 

: Myrzakulov and Sebastian! 

(20141): Elizalde et al. 1 

(2014). Investigations dealing 


with the coupling between the dark en ergy and dark 


matte r components are considered by I Brevik et al. 


( 2014 ); Timoshkin ( 20091) . A very extensive review on 
cosmological evolution with interaction between dark 
energy and dark mat ter in general is recently given by 


Bolotin et al. (|2015[) . 


A main motivation for the present study is to inves¬ 
tigate the influence from heat dissipation. We will thus 
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study the Little Rip, the Pseudo Rip, and the bounce 
cosmology, by considering a homogeneous dark fluid 
interacting with dark matter. This means that each 
fluid component constitutes a non-closed physical sys¬ 
tem, implying exchange of energy with its companion 
fluid component. Mathematically, it means the appear¬ 
ance of a source term on the right hand side of the 
energy conservation equation; cf. Eq. (HD below. We 
will assume that for each fluid component (dark energy 
and dark matter) there exists a linear and homogeneous 
equation of state; this being the simplest alternative. 
We moreover assume a spatially flat, homogeneous and 
isotropic FRW space. 

Analytical expressions will be worked out for the cos¬ 
mological cases mentioned above, and we will make ap¬ 
propriate choices for the various thermodynamical pa¬ 
rameters. Energy dissipation is conveniently expressed 
in terms of the bulk viscosity, as this kind of viscos¬ 
ity is compatible with the assumption about spatial 
isotropy. It should be mentioned, however, that this 
is actually a delicate point. The shear viscosity cal¬ 
culated in co smology on the basis o f the Boltzmann 
equation (cf. Hogeveen et al. ( 1986ll l is much larger 
than the bulk viscosity. Thus even a minute anisotropy 
in the universe might make the shear viscosity an im¬ 
portant component. The shear viscosity will however 
be omitted here, in accordance with common usage. 

Our exposition will be related to the so-called en- 
tropic force model - cf. Eq. (SD below - in which the 
equation of state has quite a general form, expressed in 
terms of powers of the Hubble parameter. 

Finally, as a general point one may ask: what are the 
physical reasons for introducing viscosity coefficients 
in the cosmic fluid? This question can be considered 
from different angles. First, from a hydrodynamical 
point of view it is apparent that the ideal (nonviscous) 
fluid model applies only under restricted circumstances. 
Once there are irreversible process involved, the intro¬ 
duction of one or two viscosity coefficients becomes in¬ 
evitable. Second, in a cosmological context the inclu¬ 
sion of viscosity broadens the applicability of the theory 
considerably. Not least so, this applies in connection 
with the transition to turbulent motion which is quite 
pro bable when the fluid a pproaches the future singular- 
(120111) ]. Even in the very early uni- 


Brevik et al. 


ity 

verse, in connection with the inflationary period, viscos¬ 
ity has been proposed in order to dea l with t he grateful 
exit; cf. iMvrzakulov and Sebastiani (l2014all . It is thus 
quite understandable that a large number of papers on 
viscosity has appeared on viscous cosmology. Readers 
interested in review papers may consu lt, for instance, 
GrOn ( 1990l) : Brevik and GrOn ( 2013 ). 


2 Dissipative cosmologies coupled with dark 
matter 

We consider a phenomenological model in which there 
occurs a constant entropic pressure associated with the 
dark matter fluid component. This falls into line with 
the so-called entropic cosmolo g|v model for a dissipa¬ 
tive u niverse (see the paper by 


Komatsu and Kimura 


(|2014l) and further references therein). This is a gener¬ 
alized form of cosmology combining bulk viscosity the¬ 
ory with creation theory of cold dark matter. 

The modified continuity and acceleration equations 
for the interacting fluid system can be written as 


p + SH{p + p) = -Q, ^ 

Pm T T Pe) — Qi 

H = -^{p+P + Pm+Pe)- 


( 1 ) 


Here H = a/a is the Hubble parameter, and = 8ttG 
with Newton’s gravitational constant G. Moreover p,p 
and Pm,Pm are the energy density and pressure respec¬ 
tively of dark energy and dark matter, and Q describes 
the interaction. A dot means derivative with respect to 
cosmic time t. We put c = 1. Quantities referring to 
present time t = 0 are given a subscript zero. 

The first l ine in Eg. (fTI) i s the co ntinuity equation for 


dark energy [Noiiri et al. (|2005aD ]. The second line is 


the continuity equation for dark matter, in which pe is 
the mentioned constant pressure term. This term de¬ 
scribes the influence from dissipation, and is connected 
with the temperature and the entropy of the Hubble 
horizon. It describes irreversible entropy, generated ei¬ 
ther by a bulk viscosity in one formulation, or with the 
creation of cold dark matter in the alternative formula- 
tion (a closer discussion on t his point can be found in 
Komatsu and Kimura (120141) and earlier papers). We 
assume that the dark matter is dust, so that its me¬ 
chanical pressure is Pm = 0 . 

The constant effective pressure Pe has to be a nega¬ 
tive quantity. Gompare here with the behavior of vis¬ 
cous cosmology where the thermodynamic pressure p is 
reduced by an amount Sf/iL, where C is the bulk viscos¬ 
ity, necessarily positive from thermodynamic grounds. 
We shall write it as Pe oc — 2iJg/A:^, with a positive con¬ 
stant of proportionality. Calling this constant {p, — 7 ) 
(because of the use of 7 in the subsequent Friedmann 
equation ([HD), we must have p > 7 > 0 . Thus 


Pe = 


2^0 . ^ 


( 2 ) 


and the Friedmann equation is written as 


H^ = -{p + Pm)+lHl 


( 3 ) 
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The additional term corresponds to an entropic- 
force term derived from revers ible entropy in the en- 


tropic cosmology model; cf. iKomatsu and Kimura 
( 20141) . The term can be interpreted as an effective 
dark energy term. The constant effective pressure Pe in 
Eq. @ consists of contributions from the entropy term 
to the pressure. They turn up in the coefficients Wi{t) 
given in Eq. (U) below. 

We consider the following general equation of state 


We now start from the gravitational equation of mo¬ 
tion 


2\H 


+ P + P + Pra+Pe = 0. 


Choosing for the parameter 7 the value 


7 = 


Pe 


A 

h'o 


( 10 ) 


( 11 ) 


for the dark energy Capozziello et al. ( 20051 ): Noiiri and and Odintsov 1 „ , ,,, . 

W ahcf Choosing for the parameter wi[t) m the equation 

of state (ID) the value 


p = -I- + W 2 {t)H'^ -I- wi{t)H, (4) 

where the thermodynamic parameters wi,W 2 ,'W 3 ,W 4 
are so far undetermined functions of the cosmic time 
t. Higher order terms are neglected because inflation in 
the early universe is not considered. 

We turn now to dissipative properties of the inter¬ 
acting dark energy universe, via proper adjustments of 
the parameters in Eq. O- 


wi{t) = - 


A 

F’ 


we can rewrite Eq. OT in the simple form 
3 


W 4 H -I- W 3 H -I- 1U2 + 




= 0 . 


( 12 ) 


(13) 


Thus, if we require that the parameters W 2 ,W 3 ,W 4 are 
related via the condition 


2.1 Little Rip model 

Let us assume the Little Rip model where the Hubble 
parameter increases exponentially with t, 

H = Hoe^*, Ho>0,X>0 (5) 

(recall that t = 0 is present time). 

Consider again the continuity equation for dark mat¬ 
ter 

Pm+^H{pm+Pe) = Q, (6) 

and insert the following expression for the interaction 
term 

Q = jHl (7) 

In this case the solution of Eq. m for the dark matter 
is given by 


Aw4 



= 0 , 


(14) 


we can from Eq. m determine the thermodynamic 
parameters as 


r W 3 {t) =-2HoW4{t)exp{Xt), 

< W 2 {t) = W 4 {t)HQ exp{2Xt) — 3/k‘^, (15) 

\ W4(t) € R. 

Insertion of these values in the equation of state o 
shows that the dark fluid pressure p acquires an expo¬ 
nentially increasing behavior when t —>■ oo. Thus the 
values for the re’s chosen in Eqs. uni) and (13 lead to 
the characteristic features of the Little Rip cosmology. 
This circumstance motivates our parameter choices, in¬ 
cluding the ansatz (HU for 7 . Physically, we have as a 
checking point that Q —t 0 when t —>■ 00 . 


2.2 Pseudo Rip model 


Pmit) = (po - Aprn) CXp ( - 


3(H - Ho) 
X 


+ Apn 


( 8 ) 


where Pm( 0 ) = po and the constant correction term 
^Pm is 


fr5 

Apm=^ 


Pe- 


( 9 ) 


Now we will investigate a cosmological model where 
the dark energy is monotonically increasing. In this 
case the Hubble parameter asymptotically approaches 
a cosmological constant, which means that the universe 
approaches a de Sitter form. 

We suppose that the Hubble parameter is given by 


Erampton et al. (l20I2fl ] 


From Eq. (| 8 ]) one sees that pm{t —>■ 00 ) = Apm- This 
correction is caused by the exchange of energy between 
dark fluid and dark matter and by dissipative properties 
in the universe. 


H = Ho- (16) 

where Ho, Hi and A are positive constants, Ho > Hi > 

0, t > 0. 
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Near the present time {t —>■ 0) we have H —>■ Hq — 
Hi, and in the late-time universe (t —>■ c») the Hubble 
parameter tends to a constant. 

We take now the coupling term to have the form 

Q = dH{pm+Pe): (17) 

where <5 is a positive nondimensional constant. This 
choice is motivated by the physically natural assump¬ 
tion that the coupling increases when the relative rate 
of change of the scale factor increases, i.e. when H 
increases. 

From Eq. (| 6 ]) we find the energy density of dark mat¬ 
ter 


Pm{t) 


{po+Pe) exp 


(<5-3) 



H-Ho +Hi 
A 


-Pe 


(18) 


with Pm{Q) = Po as before. 

The gravitational equation of motion for dark energy 
becomes 

^(i7o -H)+p+^H^+p,-^H^ = 0, (19) 

with A = ^jHq, 7 being given by Eq. (ED. 

Thus, if we choose the thermodynamic parameters 
in the equation of state (| 3 ]) to have the form 

wi = ^Ho, 

W2 = 

W3{t) = [iJiexp {-\^Hot) 

Wi{t) G R, 

we obtain the characteristics of the Pseudo Rip dissi¬ 
pative universe. Assuming that 1^4 stays finite when 
t —>■ 00 , all the other w's stay finite, and prevent the 
pressure in Eq. (jT]) from blowing up in the far future. 
Again, our choices of parameters are justified by the 
adjustment to one of the future singular cosmological 
models. 


- Ho] W4{t), 


( 20 ) 


P being a positive constant.The bounce occurs at t = 
0 (a = 1). We assume that f > 0 in the following. Now 
take the Hubble parameter to have the form 

H = 2pt, (22) 

and take the interaction between dark energy and dark 
matter to be 


Q = Qo exp(-3,5t^), (23) 

with Qo = Q(0). Thus near the bounce Q —^ Qo, 
whereas for increasing values of t, Q —)► 0 and the cou¬ 
pling becomes negligible. 

Now the solution of Eq. (jlD for the dark matter be- 

I mes 

i(.t) = {Qot + Po -kPe)exp(-3/3t^) -pe, (24) 


with pm(0) = Po as before. Just after the bounce the 
dark matter energy density thus increases linearly with 
t, whereas when t —>■ 00 , pm —>■ —Pe- The correction 
—Pe takes into account the dissipative processes. 

We make the following choice for the parameter /3, 

/3 = llHl (25) 

and write the gravitational equation for dark energy as 

^ + ( 26 ) 

The conditions ® and (El) lead to the simple form 

dMD: 


WiH'^ + woH^ -I- 



H^ + WlH +Pe 


0. (27) 


Erom Eq. (l26l) one can find the following solution 


r = 

J W2=-^, 

I woit) = -^'yH^W4{t)t, 
{. Wji(t) G R. 


(28) 


2.3 Bounce cosmology 


In the matter bounce cosmology model, the universe 
goes from an era of accelerated collapse into an era 
of expansion via a bounce acting as a singularity. Af¬ 
ter the bounce, the universe is in a matter dominated 
phase. 

Let us consider such a model where the scale factor 


a has the exponential form Bamba (|2ni3r i] 
a = exp(/3t^), 


( 21 ) 


It is seen that the singularity property at the bounce 
at t = 0 is taken care of by the coefficient wi. 

We have thus in this subsection presented a descrip¬ 
tion of dissipative bounce cosmology, related to the 
equation of state ®. 


3 Conclusion 

We have investigated some examples of Little Rip, 
Pseudo Rip, and the bounce cosmology, under condi¬ 
tions where a dark energy fluid is coupled with dark 
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matter via an interaction term Q; cf. Eq. ([T]). The 
underlying geometry is a flat Friedmann-Robertson- 
Walker universe. The energy dissipation is most natu¬ 
rally described by means of a bulk viscosity, although 
there is a mathematically equivalent formalism imply¬ 
ing creation of cold dark matter. 

The formalism that we have m ade use of, falls within 
the so -called entropic force model Komatsu and Kimura ^ 
pOldil ]. We have assumed initially that the equation 
of state for the dark fluid has the general form 0, 
in which there occurs four thermodynamic coefficients 
Wi, i G [1,4]. The various categories of singularity- 
containing universes are then obtained by appropriate 
choices of these coefficients. Corrections to the dark 
matter energy density caused by dissipative processes 
are obtained. 
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Finally, we wish to make some additional comments 
on the physical meaning of the constant pressure Pe 
introduced in Eq. o, from the viewpoint of viscous 
cosmology. As mentioned above, the modification of 
thermodynamic pressure in a viscous fluid consists in 
replacing the thermodynamic pressure p with p — 
with C the bulk viscosity. Thus, the quantity Pe cor¬ 
responds to —3H(, and the assumed constancy of Pe 
implies that has to be inversely proportional to H. 

It is worth noticing that when this condition is 
satisfied, the Friedmann formalism for the viscous 
fluid becomes easily manageable. Let us recall from 


Brevik et al. (|201l[ l how the formalism looks if one 


starts from the equation of state in the form 


p = -p- Ay/p + pe 


(29) 


for the dark fluid. Here A is a constant, with dimen¬ 
sion cm“^ in geometric units. In this case we find the 
following explicit formula for how p depends on cosmic 
time t (> 0), 


2 I Ayfp-pe 

y/ikA Ayf^-pe' 


(30) 


This can be inverted to give p as a function of t. 


P{t) = 


Pe\ I Vs 

a) ""p ^ 


kAt\+iV 


(31) 


It is seen that the universe needs an infinite time to 
reach the state p —>■ oo. This is just the characteristic 
property of Little Rip cosmology. 
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